Novel pore size-controlled, susceptibility matched, 3D-printed MRI phantoms.

PURPOSE: We report the design concept and fabrication of MRI phantoms, containing blocks of aligned microcapillaires that can be stacked into larger arrays to construct diameter distribution phantoms or fractured, to create a "powder-averaged" emulsion of randomly oriented blocks for vetting or calibrating advanced MRI methods, that is, diffusion tensor imaging, AxCaliber MRI, MAP-MRI, and multiple pulsed field gradient or double diffusion-encoded microstructure imaging methods. The goal was to create a susceptibility-matched microscopically anisotropic but macroscopically isotropic phantom with a ground truth diameter that could be used to vet advanced diffusion methods for diameter determination in fibrous tissues. METHODS: Two-photon polymerization, a novel three-dimensional printing method is used to fabricate blocks of capillaries. Double diffusion encoding methods were employed and analyzed to estimate the expected MRI diameter. RESULTS: Susceptibility-matched microcapillary blocks or modules that can be assembled into large-scale MRI phantoms have been fabricated and measured using advanced diffusion methods, resulting in microscopic anisotropy and random orientation. CONCLUSION: This phantom can vet and calibrate various advanced MRI methods and multiple pulsed field gradient or diffusion-encoded microstructure imaging methods. We demonstrated that two double diffusion encoding methods underestimated the ground truth diameter.

Limitations and generalizations of the first order kinetics reaction expression for modeling diffusion-driven exchange: Implications on NMR exchange measurements.

The study and modeling of water exchange in complex media using different applications of diffusion and relaxation magnetic resonance (MR) have been of interest in recent years. Most models attempt to describe this process using a first order kinetics expression, which is appropriate to describe chemical exchange; however, it may not be suitable to describe diffusion-driven exchange since it has no direct relationship to diffusion dynamics of water molecules. In this paper, these limitations are addressed through a more general exchange expression that does consider such important properties. This exchange fraction expression features a multi-exponential recovery at short times and a mono-exponential decay at long times, both of which are not captured by the first order kinetics expression. Furthermore, simplified exchange expressions containing partial information of the analyzed system's diffusion and relaxation processes and geometry are proposed, which can potentially be employed in already established estimation protocols. Finally, exchange fractions estimated from simulated MR data and derived here were compared, showing qualitative similarities but quantitative differences, suggesting that the features of the derived exchange fraction in this paper can be partially recovered by employing an existing estimation framework.

MRI with generalized diffusion encoding reveals damaged white matter in patients previously hospitalized for COVID-19 and with persisting symptoms at follow-up.

There is mounting evidence of the long-term effects of COVID-19 on the central nervous system, with patients experiencing diverse symptoms, often suggesting brain involvement. Conventional brain MRI of these patients shows unspecific patterns, with no clear connection of the symptomatology to brain tissue abnormalities, whereas diffusion tensor studies and volumetric analyses detect measurable changes in the brain after COVID-19. Diffusion MRI exploits the random motion of water molecules to achieve unique sensitivity to structures at the microscopic level, and new sequences employing generalized diffusion encoding provide structural information which are sensitive to intravoxel features. In this observational study, a total of 32 persons were investigated: 16 patients previously hospitalized for COVID-19 with persisting symptoms of post-COVID condition (mean age 60 years: range 41-79, all male) at 7-month follow-up and 16 matched controls, not previously hospitalized for COVID-19, with no post-COVID symptoms (mean age 58 years, range 46-69, 11 males). Standard MRI and generalized diffusion encoding MRI were employed to examine the brain white matter of the subjects. To detect possible group differences, several tissue microstructure descriptors obtainable with the employed diffusion sequence, the fractional anisotropy, mean diffusivity, axial diffusivity, radial diffusivity, microscopic anisotropy, orientational coherence (Cc) and variance in compartment's size (CMD) were analysed using the tract-based spatial statistics framework. The tract-based spatial statistics analysis showed widespread statistically significant differences (P < 0.05, corrected for multiple comparisons using the familywise error rate) in all the considered metrics in the white matter of the patients compared to the controls. Fractional anisotropy, microscopic anisotropy and Cc were lower in the patient group, while axial diffusivity, radial diffusivity, mean diffusivity and CMD were higher. Significant changes in fractional anisotropy, microscopic anisotropy and CMD affected approximately half of the analysed white matter voxels located across all brain lobes, while changes in Cc were mainly found in the occipital parts of the brain. Given the predominant alteration in microscopic anisotropy compared to Cc, the observed changes in diffusion anisotropy are mostly due to loss of local anisotropy, possibly connected to axonal damage, rather than white matter fibre coherence disruption. The increase in radial diffusivity is indicative of demyelination, while the changes in mean diffusivity and CMD are compatible with vasogenic oedema. In summary, these widespread alterations of white matter microstructure are indicative of vasogenic oedema, demyelination and axonal damage. These changes might be a contributing factor to the diversity of central nervous system symptoms that many patients experience after COVID-19.

Using deep learning to accelerate magnetic resonance measurements of molecular exchange.

Real-time monitoring and quantitative measurement of molecular exchange between different microdomains are useful to characterize the local dynamics in porous media and biomedical applications of magnetic resonance. Diffusion exchange spectroscopy (DEXSY) is a noninvasive technique for such measurements. However, its application is largely limited by the involved long acquisition time and complex parameter estimation. In this study, we introduce a physics-guided deep neural network that accelerates DEXSY acquisition in a data-driven manner. The proposed method combines sampling pattern optimization and physical parameter estimation into a unified framework. Comprehensive simulations and experiments based on a two-site exchange system are conducted to demonstrate this new sampling optimization method in terms of accuracy, repeatability, and efficiency. This general framework can be adapted for other molecular exchange magnetic resonance measurements.

Diffusion within pores fully revealed by magnetic resonance.

The diffusion propagator fully characterizes the diffusion process, which is highly sensitive to the confining boundaries and the structure within enclosed pores. While magnetic resonance has extensively been used to observe various features of the diffusion process, its full characterization has been elusive. Here, we address this challenge by employing a special sequence of magnetic field gradient pulses for measuring the diffusion propagator, which allows for "listening to the drum," mapping structural dispersity, and determining not only the pore's shape but also diffusive dynamics within it.

Cumulant expansion with localization: A new representation of the diffusion MRI signal.

Diffusion MR is sensitive to the microstructural features of a sample. Fine-scale characteristics can be probed by employing strong diffusion gradients while the low b-value regime is determined by the cumulants of the distribution of particle displacements. A signal representation based on the cumulants, however, suffers from a finite convergence radius and cannot represent the 'localization regime' characterized by a stretched exponential decay that emerges at large gradient strengths. Here, we propose a new representation for the diffusion MR signal. Our method provides not only a robust estimate of the first three cumulants but also a meaningful extrapolation of the entire signal decay.

The direction-dependence of apparent water exchange rate in human white matter.

Transmembrane water exchange is a potential biomarker in the diagnosis and understanding of cancers, brain disorders, and other diseases. Filter-exchange imaging (FEXI), a special case of diffusion exchange spectroscopy adapted for clinical applications, has the potential to reveal different physiological water exchange processes. However, it is still controversial whether modulating the diffusion encoding gradient direction can affect the apparent exchange rate (AXR) measurements of FEXI in white matter (WM) where water diffusion shows strong anisotropy. In this study, we explored the diffusion-encoding direction dependence of FEXI in human brain white matter by performing FEXI with 20 diffusion-encoding directions on a clinical 3T scanner in-vivo. The results show that the AXR values measured when the gradients are perpendicular to the fiber orientation (0.77 ± 0.13 s - 1, mean ± standard deviation of all the subjects) are significantly larger than the AXR estimates when the gradients are parallel to the fiber orientation (0.33 ± 0.14 s - 1, p < 0.001) in WM voxels with coherently-orientated fibers. In addition, no significant correlation is found between AXRs measured along these two directions, indicating that they are measuring different water exchange processes. What's more, only the perpendicular AXR rather than the parallel AXR shows dependence on axonal diameter, indicating that the perpendicular AXR might reflect transmembrane water exchange between intra-axonal and extra-cellular spaces. Further finite difference (FD) simulations having three water compartments (intra-axonal, intra-glial, and extra-cellular spaces) to mimic WM micro-environments also suggest that the perpendicular AXR is more sensitive to the axonal water transmembrane exchange than parallel AXR. Taken together, our results show that AXR measured along different directions could be utilized to probe different water exchange processes in WM.

Computing the orientational-average of diffusion-weighted MRI signals: a comparison of different techniques.

Numerous applications in diffusion MRI involve computing the orientationally-averaged diffusion-weighted signal. Most approaches implicitly assume, for a given b-value, that the gradient sampling vectors are uniformly distributed on a sphere (or 'shell'), computing the orientationally-averaged signal through simple arithmetic averaging. One challenge with this approach is that not all acquisition schemes have gradient sampling vectors distributed over perfect spheres. To ameliorate this challenge, alternative averaging methods include: weighted signal averaging; spherical harmonic representation of the signal in each shell; and using Mean Apparent Propagator MRI (MAP-MRI) to derive a three-dimensional signal representation and estimate its 'isotropic part'. Here, these different methods are simulated and compared under different signal-to-noise (SNR) realizations. With sufficiently dense sampling points (61 orientations per shell), and isotropically-distributed sampling vectors, all averaging methods give comparable results, (MAP-MRI-based estimates give slightly higher accuracy, albeit with slightly elevated bias as b-value increases). As the SNR and number of data points per shell are reduced, MAP-MRI-based approaches give significantly higher accuracy compared with the other methods. We also apply these approaches to in vivo data where the results are broadly consistent with our simulations. A statistical analysis of the simulated data shows that the orientationally-averaged signals at each b-value are largely Gaussian distributed.

On the generalizability of diffusion MRI signal representations across acquisition parameters, sequences and tissue types: Chronicles of the MEMENTO challenge.

Diffusion MRI (dMRI) has become an invaluable tool to assess the microstructural organization of brain tissue. Depending on the specific acquisition settings, the dMRI signal encodes specific properties of the underlying diffusion process. In the last two decades, several signal representations have been proposed to fit the dMRI signal and decode such properties. Most methods, however, are tested and developed on a limited amount of data, and their applicability to other acquisition schemes remains unknown. With this work, we aimed to shed light on the generalizability of existing dMRI signal representations to different diffusion encoding parameters and brain tissue types. To this end, we organized a community challenge - named MEMENTO, making available the same datasets for fair comparisons across algorithms and techniques. We considered two state-of-the-art diffusion datasets, including single-diffusion-encoding (SDE) spin-echo data from a human brain with over 3820 unique diffusion weightings (the MASSIVE dataset), and double (oscillating) diffusion encoding data (DDE/DODE) of a mouse brain including over 2520 unique data points. A subset of the data sampled in 5 different voxels was openly distributed, and the challenge participants were asked to predict the remaining part of the data. After one year, eight participant teams submitted a total of 80 signal fits. For each submission, we evaluated the mean squared error, the variance of the prediction error and the Bayesian information criteria. The received submissions predicted either multi-shell SDE data (37%) or DODE data (22%), followed by cartesian SDE data (19%) and DDE (18%). Most submissions predicted the signals measured with SDE remarkably well, with the exception of low and very strong diffusion weightings. The prediction of DDE and DODE data seemed more challenging, likely because none of the submissions explicitly accounted for diffusion time and frequency. Next to the choice of the model, decisions on fit procedure and hyperparameters play a major role in the prediction performance, highlighting the importance of optimizing and reporting such choices. This work is a community effort to highlight strength and limitations of the field at representing dMRI acquired with trending encoding schemes, gaining insights into how different models generalize to different tissue types and fiber configurations over a large range of diffusion encodings.

Q-space trajectory imaging with positivity constraints (QTI+).

Q-space trajectory imaging (QTI) enables the estimation of useful scalar measures indicative of the local tissue structure. This is accomplished by employing generalized gradient waveforms for diffusion sensitization alongside a diffusion tensor distribution (DTD) model. The first two moments of the underlying DTD are made available by acquisitions at low diffusion sensitivity (b-values). Here, we show that three independent conditions have to be fulfilled by the mean and covariance tensors associated with distributions of symmetric positive semidefinite tensors. We introduce an estimation framework utilizing semi-definite programming (SDP) to guarantee that these conditions are met. Applying the framework on simulated signal profiles for diffusion tensors distributed according to non-central Wishart distributions demonstrates the improved noise resilience of QTI+ over the commonly employed estimation methods. Our findings on a human brain data set also reveal pronounced improvements, especially so for acquisition protocols featuring few number of volumes. Our method's robustness to noise is expected to not only improve the accuracy of the estimates, but also enable a meaningful interpretation of contrast in the derived scalar maps. The technique's performance on shorter acquisitions could make it feasible in routine clinical practice.

Variance Measures for Symmetric Positive (Semi-) Definite Tensors in Two Dimensions.

Calculating the variance of a family of tensors, each represented by a symmetric positive semi-definite second order tensor/matrix, involves the formation of a fourth order tensor Rabcd. To form this tensor, the tensor product of each second order tensor with itself is formed, and these products are then summed, giving the tensor Rabcd the same symmetry properties as the elasticity tensor in continuum mechanics. This tensor has been studied with respect to many properties: representations, invariants, decomposition, the equivalence problem et cetera. In this paper we focus on the two-dimensional case where we give a set of invariants which ensures equivalence of two such fourth order tensors Rabcd and R~abcd. In terms of components, such an equivalence means that components Rijkl of the first tensor will transform into the components R~ijkl of the second tensor for some change of the coordinate system.

The sensitivity of diffusion MRI to microstructural properties and experimental factors.

Diffusion MRI is a non-invasive technique to study brain microstructure. Differences in the microstructural properties of tissue, including size and anisotropy, can be represented in the signal if the appropriate method of acquisition is used. However, to depict the underlying properties, special care must be taken when designing the acquisition protocol as any changes in the procedure might impact on quantitative measurements. This work reviews state-of-the-art methods for studying brain microstructure using diffusion MRI and their sensitivity to microstructural differences and various experimental factors. Microstructural properties of the tissue at a micrometer scale can be linked to the diffusion signal at a millimeter-scale using modeling. In this paper, we first give an introduction to diffusion MRI and different encoding schemes. Then, signal representation-based methods and multi-compartment models are explained briefly. The sensitivity of the diffusion MRI signal to the microstructural components and the effects of curvedness of axonal trajectories on the diffusion signal are reviewed. Factors that impact on the quality (accuracy and precision) of derived metrics are then reviewed, including the impact of random noise, and variations in the acquisition parameters (i.e., number of sampled signals, b-value and number of acquisition shells). Finally, yet importantly, typical approaches to deal with experimental factors are depicted, including unbiased measures and harmonization. We conclude the review with some future directions and recommendations on this topic.

B-Tensor: Brain Connectome Tensor Factorization for Alzheimer's Disease.

AD is the highly severe part of the dementia spectrum and impairs cognitive abilities of individuals, bringing economic, societal and psychological burdens beyond the diseased. A promising approach in AD research is the analysis of structural and functional brain connectomes, i.e., sNETs and fNETs, respectively. We propose to use tensor representation (B-tensor) of uni-modal and multi-modal brain connectomes to define a low-dimensional space via tensor factorization. We show on a cohort of 47 subjects, spanning the spectrum of dementia, that diagnosis with an accuracy of 77% to 100% is achievable in a 5D connectome space using different structural and functional connectome constructions in a uni-modal and multi-modal fashion. We further show that multi-modal tensor factorization improves the results suggesting complementary information in structure and function. A neurological assessment of the connectivity patterns identified largely agrees with prior knowledge, yet also suggests new associations that may play a role in the disease progress.

Enforcing necessary non-negativity constraints for common diffusion MRI models using sum of squares programming.

In this work we investigate the use of sum of squares constraints for various diffusion-weighted MRI models, with a goal of enforcing strict, global non-negativity of the diffusion propagator. We formulate such constraints for the mean apparent propagator model and for spherical deconvolution, guaranteeing strict non-negativity of the corresponding diffusion propagators. For the cumulant expansion similar constraints cannot exist, and we instead derive a set of auxiliary constraints that are necessary but not sufficient to guarantee non-negativity. These constraints can all be verified and enforced at reasonable computational costs using semidefinite programming. By verifying our constraints on standard reconstructions of the different models, we show that currently used weak constraints are largely ineffective at ensuring non-negativity. We further show that if strict non-negativity is not enforced then estimated model parameters may suffer from significant errors, leading to serious inaccuracies in important derived quantities such as the main fiber orientations, mean kurtosis, etc. Finally, our experiments confirm that the observed constraint violations are mostly due to measurement noise, which is difficult to mitigate and suggests that properly constrained optimization should currently be considered the norm in many cases.

Using the Wild Bootstrap to Quantify Uncertainty in Mean Apparent Propagator MRI.

Purpose: Estimation of uncertainty of MAP-MRI metrics is an important topic, for several reasons. Bootstrap derived uncertainty, such as the standard deviation, provides valuable information, and can be incorporated in MAP-MRI studies to provide more extensive insight. Methods: In this paper, the uncertainty of different MAP-MRI metrics was quantified by estimating the empirical distributions using the wild bootstrap. We applied the wild bootstrap to both phantom data and human brain data, and obtain empirical distributions for the MAP-MRI metrics return-to-origin probability (RTOP), non-Gaussianity (NG), and propagator anisotropy (PA). Results: We demonstrated the impact of diffusion acquisition scheme (number of shells and number of measurements per shell) on the uncertainty of MAP-MRI metrics. We demonstrated how the uncertainty of these metrics can be used to improve group analyses, and to compare different preprocessing pipelines. We demonstrated that with uncertainty considered, the results for a group analysis can be different. Conclusion: Bootstrap derived uncertain measures provide additional information to the MAP-MRI derived metrics, and should be incorporated in ongoing and future MAP-MRI studies to provide more extensive insight.

Orientationally-averaged diffusion-attenuated magnetic resonance signal for locally-anisotropic diffusion.

Diffusion-attenuated MR signal for heterogeneous media has been represented as a sum of signals from anisotropic Gaussian sub-domains to the extent that this approximation is permissible. Any effect of macroscopic (global or ensemble) anisotropy in the signal can be removed by averaging the signal values obtained by differently oriented experimental schemes. The resulting average signal is identical to what one would get if the micro-domains are isotropically (e.g., randomly) distributed with respect to orientation, which is the case for "powdered" specimens. We provide exact expressions for the orientationally-averaged signal obtained via general gradient waveforms when the microdomains are characterized by a general diffusion tensor possibly featuring three distinct eigenvalues. This extends earlier results which covered only axisymmetric diffusion as well as measurement tensors. Our results are expected to be useful in not only multidimensional diffusion MR but also solid-state NMR spectroscopy due to the mathematical similarities in the two fields.

Multimodal integration of diffusion MRI for better characterization of tissue biology.

The contrast in diffusion-weighted MR images is due to variations of diffusion properties within the examined specimen. Certain microstructural information on the underlying tissues can be inferred through quantitative analyses of the diffusion-sensitized MR signals. In the first part of the paper, we review two types of approach for characterizing diffusion MRI signals: Bloch's equations with diffusion terms, and statistical descriptions. Specifically, we discuss expansions in terms of cumulants and orthogonal basis functions, the confinement tensor formalism and tensor distribution models. Further insights into the tissue properties may be obtained by integrating diffusion MRI with other techniques, which is the subject of the second part of the paper. We review examples involving magnetic susceptibility, structural tensors, internal field gradients, transverse relaxation and functional MRI. Integrating information provided by other imaging modalities (MR based or otherwise) could be a key to improve our understanding of how diffusion MRI relates to physiology and biology.

Bayesian uncertainty quantification in linear models for diffusion MRI.

Diffusion MRI (dMRI) is a valuable tool in the assessment of tissue microstructure. By fitting a model to the dMRI signal it is possible to derive various quantitative features. Several of the most popular dMRI signal models are expansions in an appropriately chosen basis, where the coefficients are determined using some variation of least-squares. However, such approaches lack any notion of uncertainty, which could be valuable in e.g. group analyses. In this work, we use a probabilistic interpretation of linear least-squares methods to recast popular dMRI models as Bayesian ones. This makes it possible to quantify the uncertainty of any derived quantity. In particular, for quantities that are affine functions of the coefficients, the posterior distribution can be expressed in closed-form. We simulated measurements from single- and double-tensor models where the correct values of several quantities are known, to validate that the theoretically derived quantiles agree with those observed empirically. We included results from residual bootstrap for comparison and found good agreement. The validation employed several different models: Diffusion Tensor Imaging (DTI), Mean Apparent Propagator MRI (MAP-MRI) and Constrained Spherical Deconvolution (CSD). We also used in vivo data to visualize maps of quantitative features and corresponding uncertainties, and to show how our approach can be used in a group analysis to downweight subjects with high uncertainty. In summary, we convert successful linear models for dMRI signal estimation to probabilistic models, capable of accurate uncertainty quantification.

Influence of the size and curvedness of neural projections on the orientationally averaged diffusion MR signal.

Neuronal and glial projections can be envisioned to be tubes of infinitesimal diameter as far as diffusion magnetic resonance (MR) measurements via clinical scanners are concerned. Recent experimental studies indicate that the decay of the orientationally-averaged signal in white-matter may be characterized by the power-law, E(q) ∝ q-1, where q is the wavenumber determined by the parameters of the pulsed field gradient measurements. One particular study by McKinnon et al. [1] reports a distinctively faster decay in gray-matter. Here, we assess the role of the size and curvature of the neurites and glial arborizations in these experimental findings. To this end, we studied the signal decay for diffusion along general curves at all three temporal regimes of the traditional pulsed field gradient measurements. We show that for curvy projections, employment of longer pulse durations leads to a disappearance of the q-1 decay, while such decay is robust when narrow gradient pulses are used. Thus, in clinical acquisitions, the lack of such a decay for a fibrous specimen can be seen as indicative of fibers that are curved. We note that the above discussion is valid for an intermediate range of q-values as the true asymptotic behavior of the signal decay is E(q) ∝ q-4 for narrow pulses (through Debye-Porod law) or steeper for longer pulses. This study is expected to provide insights for interpreting the diffusion-weighted images of the central nervous system and aid in the design of acquisition strategies.

Asymmetric Orientation Distribution Functions (AODFs) revealing intravoxel geometry in diffusion MRI.

Characterization of anisotropy via diffusion MRI reveals fiber crossings in a substantial portion of voxels within the white-matter (WM) regions of the human brain. A considerable number of such voxels could exhibit asymmetric features such as bends and junctions. However, widely employed reconstruction methods yield symmetric Orientation Distribution Functions (ODFs) even when the underlying geometry is asymmetric. In this paper, we employ inter-voxel directional filtering approaches through a cone model to reveal more information regarding the cytoarchitectural organization within the voxel. The cone model facilitates a sharpening of the ODFs in some directions while suppressing peaks in other directions, thus yielding an Asymmetric ODF (AODF) field. We also show that a scalar measure of AODF asymmetry can be employed to obtain new contrast within the human brain. The feasibility of the technique is demonstrated on in vivo data obtained from the MGH-USC Human Connectome Project (HCP) and Parkinson's Progression Markers Initiative (PPMI) Project database. Characterizing asymmetry in neural tissue cytoarchitecture could be important for localizing and quantitatively assessing specific neuronal pathways.